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ABSTRACT 

We report the detection of the Hi hne at 21 cm in the direction of a Ori with the 
Nangay Radiotelescope and with the Very Large Array. 

The observations confirm the previous detection of Hi emission centered on a Ori, 
but additionally reveal for the first time a quasi-stationary detached shell of neutral 
atomic hydrogen ^ 4' in diameter (0.24 pc at a distance of 200 pc). The detached 
shell appears elongated in a direction opposite to the star's space motion. 

A simple model shows that this detached atomic gas shell can result from the 
collision of the stellar wind from a Ori with the local interstellar medium (ISM) . It 
implies that a Ori has been losing matter at a rate of ^ 1.2 x 10^^ Mq yr^^ for the 
past 8 X 10^ years. 

In addition, we report the detection of atomic hydrogen associated with the far- 
infrared arc located 6' north-east of a Ori, that has been suggested to trace the bow 
shock resulting from the motion of the star through the ISM. We report also the 
detection by the Galaxy Evolution Explorer (GALEX) of a far-UV counterpart to this 
arc. 
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1 INTRODUCTION 

Evolved late-type stars (AGB stars and red supergiants) are 
observed to undergo mass loss at high rates. They are thus 
surrounded by expanding circumstellar envelopes (CSEs). 
Observations show a wide variety of different properties of 
the central stars and of the CSEs, which are not always easily 
related. The time variability of the mass loss phenomenon 
may explain why it is difficult to relate the properties of 
CSEs to those of the central stars. Also, some stars may 
share the same observational properties, but be in slightly 
different evolutionary stages, with different initial masses 
and different mass loss histories. It would be extremely use- 
ful to get overviews of CSEs that would allow us to smooth 
the chaotic effects of the mass loss time-variability and as 
well to measure the total mass expelled by the central stars. 
With this goal in mind, we are exploring the properties of 
CSEs around evolved stars by using the atomic hydrogen 
(Hi) line at 21 cm. 

Indeed, hydrogen is expected to be the dominant species 
in stellar winds of evolved stars. It is also expected to be in 
atomic form right from the stellar atmosphere if the effec- 
tive temperature. Tea, is larger than 2500 K, and to remain 
atomic in the expanding wind (Glassgold & Huggins 1983) . 



Also, atomic hydrogen should be protected from photoion- 
ization by the abundant atomic hydrogen in the interstellar 
medium (ISM) which absorbs UV photons of energy larger 
than 13.6 eV. The Hi line at 21 cm should thus be an im- 
portant tracer of stellar winds for late- type stars with Tcfi> 
2500 K, including the outer CSEs where molecules are dis- 
sociated. Hi studies should bring unique information on the 
kinematics in these media, and in particular in the regions 
of interaction with the ISM. 

However the detection of red giants in the H I line at 21 
cm is hampered by the competing emission from atomic hy- 
drogen in the ISM (which ironically protects the circumstel- 
lar atomic hydrogen from destruction) along the same lines 
of sight. For a long time this difficulty precluded Hi obser- 
vational studies of evolved stars. However, recently, taking 
advantage of upgraded instruments, we addressed this topic 
and obtained several detections of AGB stars (Gerard & 
LeBertre 2006, Matthews & Reid 2007). In general, a nar- 
row emission line with a quasi-gaussian profile (FWHM ~ 
3kms~^) is observed close to the stellar radial velocity (V*). 

The fate of massive stars is strongly dependent on their 
mass loss history, which is still poorly understood, in par- 
ticular when they are red supergiants (van Loon 2010). Our 
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Figure 1. Evolutionary tracks from Meynet &; Maeder (2003; 
Z = 0.02) for stars of initial masses 15 and 20 Mq until the end of 
He-burning; full lines : models with rotation (Vi = 300 kms"'^), 
dashed lines: models without rotation. The star symbol (*) marks 
the position of a Ori assuming a distance of 200 pc (log L/Lq = 
5.1; Teff= 3650 K). 



previous successes encouraged us to extend H I observations 
to this class of stars. Although suffering from strong galactic 
H I emission, the supergiant a Ori was detected by Bowers & 
Knapp (1987) using the Very Large Array. At the position of 
the star, they found emission from the front and back sides 
of a CSE expanding at a velocity of 14 kms"'^. 

In this paper, we report new 21-cm observations of a Ori 
that we analyze in combination with the previous data from 
Bowers & Knapp (1987), in order to revisit in Hi this oth- 
erwise well-documented red supergiant. 



2 BASIC DATA 

a Ori (Betelgeuse) is the prototype of red supergiants. It 
has a spectral type M2Iab. Perrin et al. (2004) have mea- 
sured its effective temperature by infrared interferometry: 
Teft=3641±53K. This determination fits well with the re- 
vised temperature scale of Galactic Red Supergiants by 
Levesque et al. (2005, 3650 K for a Ori). It indicates clearly 
that atomic hydrogen should be the dominant species in the 
atmosphere and outflow of a Ori. 

Although it is the closest supergiant, its distance is not 
well known. Hipparcos measured a parallax of 7.63±1.64 
(Ferryman et al. 1997). It has been revised to 6.55±0.83 
mas (van Leeuwen 2007), which translates to a distance of 
153 pc (135-175 pc). But Harper et al. (2008), using multi- 
wavelength radio data obtained at the Very Large Array 
(VLA), derived a smaller parallax/larger distance. When 
combining the original Hipparcos results with VLA data, 
they get d=197±45 pc. Harper et al. find also proper mo- 
tions which are reduced with respect to Hipparcos. They 
suggest that the Hipparcos parallax is in error because of 
the large angular size of the star and the motions of the stel- 
lar photocenter due to convection (Lim et al. 1998, Haubois 
et al. 2009, Ohnaka et al. 2011). In the following we follow 
Harper et al. (2008) and adopt a nominal distance of 200 pc. 
At this distance, the bolometric luminosity corresponds to 
log L/L0« 5.10 ± 0.22, which, from the Meynet & Maeder 




Figure 2. Upper panel: IRAS image at 60 fim centered on a 
Ori (star symbol). Note the arc (which is pointed by the three 
arrows), ~ 6' north-east of a Ori, and the bar, 10' north-east 
(Noriega-Crespo et al. 1997). Lower panel: the same field observed 
by GALEX. Note the faint narrow counterpart to the 60 /^m arc 
detected by IRAS from PA ~ 0° to 100°. The scale is in counts 
per second and per pixel (1.5"). 



(2003) models (see Fig. [T]), indicates a progenitor that was 
an 09V main-sequence star of 20 M0 (Harper et al. 2008). 
Due to mass loss from the main sequence to the red super- 
giant stage, the present stellar mass would be ~ 12-16 M© . 
Recently using predictions of limb darkening by stellar at- 
mosphere models and high spatial resolution observations 
at 1.64 /^m, Neilson et al. (2011) estimate a present mass 
of 11.6l3gM0, in good agreement with Maeder & Meynet's 
2OM0-models. However, this should not be overemphasized 
since, as mentioned above, the mass loss history of red su- 
pergiants is rather insecure. 

The proper motions, revised by Harper et al. (2008), 
are 24.95 masyr"^ in right ascension, and 9.56 masyr"^ 
in declination. For a distance of 200 pc, and assuming 
a solar motion, Uo = 11.1 kms~^, Vo = 12.24kms~ , and 
Wo = 7.25kms~^ (Schonrich et al. 2010), we derive a ve- 
locity in the plane of the sky of 30kms~^ with an apex at 
PA=53°. 

Like most supergiants, a Ori is surrounded by an ex- 
tended circumstellar envelope. The physical processes re- 
sponsible for red supergiant mass-loss are unclear, but may 
be related to atmospheric convection (Lim et al. 1998, 
Ohnaka et al. 2011). The convective cells could levitate cool 
gas above the photosphere and lead to mass loss (Josselin 
& Flez 2007, Chiavassa et al. 2011). KerveUa et al. (2011) 
have obtained spectacular infrared images with the Very 
Large Telescope of the European Southern Observatory that 
show an inhomogeneous spatial distribution of the material 
ejected by the central star at scales corresponding to several 
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stellar radii (0.5-3"). IRAS (Infrared Astronomical Satellite) 
detected an extended structure at 60 /im, which Young et 
al. (1993a) fitted with a detached shell model of inner ra- 
dius 0.2' (0.01 pc at 200 pc) and outer radius 10.5' (0.6 pc). 
However, through a careful analysis of high-resolution IRAS 
data, Noriega-Crespo et al. (1997) discovered a ring ('far- 
IR arc') of mean radius 6', mostly confined north-east of a 
Ori around PA= 60±10°, in good agreement with the apex 
of the proper motion (Fig. [SJ . They interpret this structure 
as a bow shock resulting from the interaction of the super- 
giant wind with the ISM. Assuming that the IR emission is 
due to dust, Noriega-Crespo et al. evaluate the mass in the 
far-IR arc to be ~ 0.034 Mq (at 200 pc). This bow-shock 
structure was confirmed by Akari at far-infrared wavelengths 
(Ueta et al. 2008). A straight filament ('bar') is seen at 10' 
from the central star ahead of the bow shock. Although per- 
pendicular to the direction of motion and of maximum in- 
tensity close to the apex of the bow shock, it does not seem 
related to a Ori in a simple manner, but may have con- 
tributed to the extended emission reported by Young et al. 
(1993a), which thus appears to be a blend of two distinct 
sources. 

Very recently, a Ori was observed by Herschel at 70 /im 
and 160 /xm (Cox et al. 2012). The images have an un- 
precedented spatial resolution (6" and 12", respectively) and 
bring many new details that could not be seen on IRAS and 
Akari images. In particular, they reveal sub-structures in 
the far-IR arc discovered by IRAS (upper panel of Fig. [2]), 
which breaks up into at least three thin shells. It is also 
clear that these shells are not circular, and that the radii of 
curvature, close to apex, are larger than the distance to a 
Ori, as would be expected from a bow shock. There are also 
hints of additional arcs at smaller radii. 

Finally, our present analysis has uncovered a faint coun- 
terpart to this ring in far-ultraviolet (FUV ~154 nm) imag- 
ing obtained by GALEX (Galaxy Evolution Explorer; Mar- 
tin et al. 2005). The GALEX data presented in Fig. [3 were 
obtained, from Nov. 29, 2008, to Jan. 14, 2009, as part of 
the Nearby Galaxies Survey (tile 5786) with an exposure 
time of 48915 sec. The GR6 pipeline calibrated data were 
retrieved from the GALEX archive (Morrissey et al. 2007), 
and smoothed by a Gaussian with kernel width 3 pixels 
(pixel = 1.5"). The FUV feature is coincident (in projection) 
with the outer shell revealed by Herschel. It can be seen most 
clearly north and northeast of the star, although there are 
hints that the entire structure subtends an angle of at least 
~105° (see arrows on Fig. [2|). In the near-UV (~ 232 nm) 
the image is unfortunately saturated by a ghost of the cen- 
tral star. A similar, narrow and arc-shaped, UV-feature has 
been reported around IRC -1-10216 by Sahai & Chronopou- 
los (2010). They ascribed the emission to coUisionally ex- 
cited H2 molecules, although scattering of the interstellar 
radiation field by small dust particles could as well be con- 
sidered. On the other hand, the a Ori bar does not seem to 
be detected by GALEX. To our knowledge the far-UV emis- 
sion associated with the far-IR arc of a Ori has not been 
reported before. 

The circumstellar envelope of a Ori has also 
been studied in the CO rotational lines. Huggins 
(1987) observed a CO(2-l) emission line centered at 
V*=3.7±0.4kms~^ (LSR). It implies a 3D space velocity of 
31 kms~^, mostly in the plane of the sky. From the profile. 
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Figure 3. NRT frequency-switclied spectrum obtained on the 
position of a Ori (continuous line) compared to the LAB spectrum 
(dashed line, Kalberla et al. 2005). We adopted a main beam 
efficiency of 0.60 for the NRT at 21 cm. The vertical line marks 
the CO centroid velocity of a Ori and the horizontal bar, the 
velocity extent of its CO emission (Huggins 1987). 



Huggins derives an expansion velocity, Vcxp=14.3kms 

and a mass-loss rate of 

pc). 



-1x10"** Mgyr"^ (scaled to 200 



3 NRT SINGLE-DISH OBSERVATIONS 

The Nangay Radiotelescope (NRT) is a meridian tele- 
scope with a rectangular aperture of effective dimensions 
160 mx 30 m. The beam-width at 21 cm is 4' in right as- 
cension and 22' in declination. The point source efficiency 
is 1.4 K Jy~^, and the system noise temperature, 35 K. The 
observations of a Ori have been obtained in the course of a 
large survey programme (Gerard et al. 2011) over the period 
from October 2006 to September 2011. 

In Fig. |3] we present a spectrum obtained in the 
frequency-switched mode of observation. It illustrates that 
in the direction of a Ori the Hi emission is dominated by 
galactic emission due to interstellar atomic hydrogen along 
the same line-of-sight. The emission peaks at ^ 7.5 kms~^. 
Our spectrum is in good agreement with the spectrum 
extracted from the Leiden/Bonn/ Argentina (LAB) survey 
(Kalberla et al. 2005; ~ 36'). No distinctive feature can 
be seen in the velocity range expected from the CO obser- 
vations of a Ori (Huggins 1987) . 

Data were also obtained in the position-switched mode 
of observation with the central star placed at the "on- 
position", and the "off-positions" taken at ±2', ±4', and 
±6' in the east-west direction (Fig. |4]). These spectra are 
strongly affected by interstellar confusion. Nevertheless, an 
emission feature is detected at the radial velocity expected 
for a Ori (3 kms~^). This is seen in the spectrum with the 
smallest throw (2'= 1/2 beam in right ascension, left panel) 
with a peak intensity of ^ 0.5 Jy. The intensity is about 
1 Jy at ±4' (1 beam, middle panel) and about 1.2 Jy at ±6' 
(1.5 beam, right panel). However, the confusion increases 
rapidly with throw (see for instance the emission feature 
at -2kms~^, and the absorption troughs at -1-10 and -|-20 
kms~^, that may affect the a Ori line-profiles). We have 
also obtained spectra with larger throws, but the confusion 
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becomes so large that the emission feature at 3 kms~^ gets 
drowned. 

FinaUy, position-switched spectra were obtained with 
the on-positions at -1-11' (1/2 beam north) and at -11' (1/2 
beam south) of the central star. We find that the H I source is 
centered on the stellar position in declination and that there 
is no evidence of an extension in the north-south direction 
(FWHM ^11'). 

Despite the high level of confusion, we can infer that 
an Hi source, coincident with q Ori, is detected at Visr= 3 
kms~^. This source does not appear extended in right ascen- 
sion as compared to the NRT beam (diameter ~ 4'). The line 
profile is narrow with a FWHM ~ 3 kms~^. The integrated 
intensity is ~ 5 Jyxkms"^, which translates to ~ 0.05 M© in 
atomic hydrogen at 200 pc. This relatively strong emission 
in Hi at a velocity for which Bowers & Knapp (1987) did not 
report any signal was puzzling and motivated a reassessment 
of the case at the Very Large Array (VLA). 



4 VLA OBSERVATIONS 

4.1 C Configuration Observations 

Bowers & Knapp (1987) previously observed a Ori in the H I 
21-cm line with the VLA of the National Radio Astronomy 
Observatory (NRAO) in its C configuration (0.08 to 3.2 km 
baselines) Ij These data comprise 32 frequency channels in 
a single (right circular) polarization with a channel spacing 
of 6.10 kHz (1.29 kms~^) after on-line Hanning smoothing, 
resulting in a total bandwidth of 0.189 MHz (~41 kms"^). 
The band was centered at an LSR velocity of -f5.0 kms^^. 
In total, ~5.5 hours of integration time was spent on-source. 
Further details on the observations can be found in Bowers 
& Knapp (1987). 

We downloaded the Bowers & Knapp data from the 
VLA archive and reprocessed them, following standard pro- 
cedures for data flagging and gain calibration within the As- 
tronomical Image Processing System (AIPS). The flux den- 
sities of our calibration sources are summarized in Table [T] 
No bandpass calibration was possible owing to contamina- 
tion of the bandpass calibrator by H i absorption across more 
than half of the band. 



4.2 D Configuration Observations 

To obtain improved sensitivity to large-scale H I emission to- 
ward a Ori, we obtained additional VLA observations of the 
star using the VLA D configuration (0.035-1.0 km baselines) 
over the course of two observing sessions on 2010 January 
2-3 (4 hours) and 2010 January 3-4 (6 hours). In total, ~6.7 
hours were spent on-source. 

The VLA correlator was used in lA mode with a 
1.56 MHz bandpass, yielding 512 spectral channels with 
3.05 kHz (~0.64 kms~^) spacing in a single (right circu- 
lar) polarization. The band was centered at an LSR velocity 
of 3.7 kms~^. Observations of a Ori were interspersed with 
observations of a phase calibrator, 0632+103, approximately 
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every 20 minutes. 3C48 (0137-1-331) was used as a flux cal- 
ibrator, and an additional strong point source (0521-1-166) 
was observed as a bandpass calibrator (Table [iJ. To insure 
that the absolute flux scale and bandpass calibration were 
not corrupted by Galactic emission in the band, the flux 
and bandpass calibrators were each observed twice, with fre- 
quency shifts of -1-1.27 MHz and —1.19 MHz, respectively, 
relative to the band center used for the observations of a Ori. 
0632+103 was also observed once at each of these offset fre- 
quencies to permit more accurate bootstrapping of the ab- 
solute flux scale to the a Ori data. We estimate that the 
resulting absolute flux scale has an uncertainty of ~10%. 

Data processing was performed using AIPS. At the 
time of our observations, the VLA contained 21 opera- 
tional antennas with L-band receivers, 19 of which had 
been retrofitted as part of the Expanded Very Large Array 
(EVLA) upgrade. Data obtained during this EVLA tran- 
sition period require special care during calibration|j After 
applying the latest available corrections to the antenna posi- 
tions and performing an initial excision of corrupted data, we 
computed and applied a bandpass calibration to our spectral 
line data to remove closure errors on VLA-EVLA baselines. 
The bandpass was normalized using channels 164-447, thus 
excluding the portion of the band affected by aliasing. We 
next computed a frequency-averaged (channel 0) data set for 
use in calibrating the frequency-independent complex gains, 
again using channels 164-447. Following gain calibration, we 
applied time-dependent frequency shifts to the data to com- 
pensate for changes caused by the Earth's motion during 
the course of the observations. At this stage, we also ap- 
plied Hanning smoothing in velocity and discarded every 
other channel, yielding a 256 channel data set with a ve- 
locity resolution of ~1.3 kms^^. There were no continuum 
sources brighter than ~0.02 Jy within our primary beam and 
none detected within a single spectral channel within ~ 10' 
of the stellar position, hence no continuum subtraction was 
performed. 



4.3 A Combined VLA Data Set 

To obtain a data set with the maximum possible u-v cover- 
age, we have concatenated the VLA C and D configuration 
data over the velocity range where the two data sets overlap 
(—14.3 ^ Vlsr ^ 24.3 kms~^). This velocity range is suf- 
ficient to encompass the range of CO velocities previously 
observed toward a Ori (see Fig. [3]). The u, v, and w values 
of the C configuration data were precessed to J2000.0 values 
prior to concatenation. 

We imaged the combined data using the AIPS IMAGR 
program with 'robust +1' weighting of the visibilities. The 
large-scale Galactic emission along the line-of-sight (Fig. [Sjl 
is poorly spatially sampled by the VLA, even in the com- 
bined C+D configuration data set, resulting in patterns of 
strong positive and negative mottling across a number of 
the channel images. Because such emission cannot be readily 
deconvolved by standard CLEAN deconvolution algorithms, 
we performed only a shallow clean (20 iterations) on each 
spectral channel to reduce sidelobes from strong compact 
emission features. The rms noise per channel is variable and 
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Figure 4. a Ori spectra obtained witii tiie NRT in tiie position-switchied mode. Left: off-positions taken at ±2'. Centre: off-positions 
taken at ±4'. Right: off-positions taken at ±6'. The thin lines correspond to the modeled spectra discussed in Sect. [5] 



Table 1. VLA Calibration Sources 



Source 



q(J2000.0) 5(J2000.0) Flux Density (Jy) Date 



3048=^ 

0521-1-166'' 

0532-1-075'^ 
0632-1-103'= 



01 37 41.29 4-33 09 35.13 



05 21 09.88 16 38 22.05 



05 32 38.99 

06 32 15.32 



07 32 43.34 
10 22 01.67 



16.038* 

15.8778* 

8.66±0.08t 

8.40±0.12t 

2.11±0.01 

2.39±0.05t 

2.25±0.04t 



1985 Sep 22 
2010 Jan 2-4 
2010 Jan 2-3 
2010 Jan 3-4 
1985 Sep 22 
2010 Jan 2-3 
2010 Jan 3-4 



Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 

* Adopted flux density at 1420.76 MHz. 
t Mean computed flux density from observations at 1419.57 and 1422.0 MHz; see §|4] 

a Primary flux calibrator 

b Bandpass calibrator 

c Phase calibrator 



in most channels is dominated by the complex, large-scale 
emission. 

Despite these complications, we find clear evidence of 
Hi emission associated with a Ori. The channel maps from 
-11.8 to -f23.0 kms~^ are presented in Fig. \5\ They are 
strongly affected by extended interstellar emission for veloc- 
ities above -5kms~'^. Nevertheless an isolated feature coin- 
cident with a Ori is clearly seen from —10.5 to — S.Skms^'^. 
In order to facilitate the identification of small sources, we 
present in Fig. [S] the channel maps obtained when selecting 
only the largest u-v spacings (^ 0.4 kA). The spatial reso- 
lution is ~ 34". A source coincident with the central star 
is detected in the channel maps from -9.2 to -6.6 kms~^, 
and from -1-17.9 to -1-19.2 kms~^. These detections are ob- 
tained in channels which correspond approximately to the 
extreme velocities expected from the wind of a Ori (V*- 



Vo 



-10.6kms-\ V*-hVo 



18.0 km s" 



We find also emission associated with a Ori at 
-1-2.4 kms~^ and -1-3.7 kms^'^. Enlarged versions of the 
two channel maps are shown in Fig. [7] (all baselines) and 
in Fig. [8] (limited to baselines larger than 0.2 kA). At -1-2.4 
kms"'^ (upper panels), the emission peaks on the central 
star. On the other hand, in the -1-3.7 kms"'^ channel map, 
one notes rather a ring of emission peaks surrounding the 
position of the star (bottom panels in Figs. [7| and [8]). Out- 
side this ring of ~ 4' in size the emission decreases steeply, 
except south-west where a plateau of emission seems to ex- 
tend away from the star. The Hi emission associated with 



Q Ori in these two channels fits well with that observed at 
the NRT (Sect. [3]). 

In addition, the channel maps from -1-6.3 to 
-1-10.1 kms~^ (Fig. [6| show Hi emission seemingly associ- 
ated with the arc of far- infrared emission discovered by IRAS 
6' north-east of a Ori (Noriega-Crespo et al. 1997). Several 
clumps are also noticeable along the far-IR arc at lower ve- 
locities, from -1.5kms~^ to +5kms~^, and at -1-10 km s~^. 
In Fig. [9] the Hi emission integrated over the range -1.5 to 
8.9 kms"'^ is overlaid on the IRAS image. The integrated 
intensity of the H I emission associated with the far-IR arc 
is ~4.9 Jyxkms^'^. This is rather uncertain since the struc- 
tures blend with background emission at low levels and be- 
cause the maps in Fig. [5] are probably missing some fiux. 
It translates to a mass of ~ 0.047 M© in atomic hydrogen. 
Assuming 10% in helium, the total mass can be estimated 
to be ~ 0.068 M0, which is within a factor 2 of the Noriega- 
Crespo et al. (1997) estimate for the mass in the far-IR arc. 
There is perhaps also Hi emission at +7.6kms~^ associated 
with the bar ahead of the arc, but the spatial coincidence is 
less convincing. 

Based on their C configuration data. Bowers & Knapp 
(1987) reported emission associated with a Ori at ~ 
— 9kms~^ and at ~-|-17kms~^, but not at ~-f3 kms~^. 
The reason is that they presented results based only on the 
largest u-v spacings (Js 0.4 kA). By doing so, they filtered 
out emission on large scales. As they noted, it is an effi- 
cient approach to get rid of most of the confusing interstel- 
lar emission. However, excluding the short spacings from the 
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Figure 5. Channel maps obtained at the VLA by combining the data in configuration C of Bowers & Knapp (1987) with our 2010 data 
in configuration D. The data have a "robust +1" weighting, and the synthesized beam has a size of 37.4"x32.0". The contour levels 
are (-11.2, -8, -5.6, -4, -2.8, -2, 2..., 11.2)x3.9 mjy/beam. The negative contours are in gray, and the positive contours in red. The star 
symbol marks the position of a Ori. 
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Figure 7. Channel maps at 2.4 km s ^ (top) and at 
3.7kms~^ (bottom). The star symbol marks the position of a 
Ori. The contour levels are (-11.2, -8, -5.6, -4, -2, 2,..., 11.2)x3.9 
mjy/beam. Dashed lines represent negative contours. Note, in the 
map at 3.7 kms~^, the emission peaks that surround the central 
star and appear to define the a Ori gas detached shell. 



C-configuration data removes also all signatures of the shell. 
In particular for a Ori at ~-|-3knis~^, it removed the source 
of ~ 4' in size associated with a Ori which was detected by 
the NRT. If the same filtering is applied to the merged C+D 
data, only the ring of emission peaks is left, as can be seen 
in the 3.7 kms~^ channel map of Fig. [S] 

Conversely, when we isolate the data obtained in con- 
figuration D, we only detect unambiguously a point source 
around -8 km s~^ . The channel maps at LSR velocities larger 
than -5kms~^ suffer heavily from interstellar emission. In 
spite of this confusion, a source is detected at +3kms^^, 
but it does not show the complex structure visible in Fig. [7] 
because of the lower spatial resolution. Also a peak of emis- 
sion is visible around -|-17kms~^, but it is immersed in a 
larger structure. The combination of data obtained in con- 
figurations C and D is thus essential to reveal the full range 
of features that have been described above. 



Figure 8. Same as Fig. Q but restricted to baselines larger than 
0.2kA. The contour levels here are (-8, -5.6, -4, 4, 5.6, 8, 11.2)x3.9 
mJy/beam. 



5 INTERPRETATION 

A star which is undergoing mass loss gets progressively 
surrounded by an expanding circumstellar shell. However, 
the collision of the supersonic stellar wind with the sur- 
rounding medium (ISM or an older envelope) produces a 
shell of denser material that is slowly expanding (Lamers & 
Cassinelli 1999). "Detached shells" have been reported to be 
detected by IRAS around many evolved stars (Young et al. 
1993a). The emission of the detached shells observed at 60 
and 100 /im is commonly assumed to be continuum emission 
radiated thermally by dust (Kerschbaum et al. 2010). Young 
et al. (1993b) developed a model of a dust shell expanding 
into the ISM that simulates the formation of a detached shell 
and accounts for the IRAS results. 

The Young et al. (1993b) scenario was revived by Lib- 
ert et al. (2007) who developed a spherical model adapted 
to the gas and to the modeling of its Hi line emission. In 
this model a detached shell is produced by a stellar outfiow 
that is abruptly slowed down at a termination shock. This 
first shock defines the inner limit of the detached shell (ri). 
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Figure 9. Hi emission from Fig. |6] summed from -1.5 to 8.9 km/s. 
The contour levels are (-8., -5.6, -4., -2.8, 2.8, 4., 5.6, 8.)xl.8 
mjy/beam. Tlie background image is from IRAS (see Fig. [2]l. 
Note the association between H I emission and the far-IR arc. 



The outer limit (r2) is defined by the leading shock (bow 
shock) where external matter is compressed by the expand- 
ing shell (Fig, nop . Between these two limits, the detached 
shell is composed of compressed circumstellar matter and 
interstellar matter separated by a contact discontinuity at 
r/ . The circumstellar matter is decelerated and heated when 
crossing the shock in ri. As it cools down from T^ (in ri) 
to T/ (in r/), the expansion velocity decreases further, from 
V]^ to v/ , and the density increases from n]*" to n7 . Adopting 
an arbitrary temperature profile, Libert et al. (2007) solve 
numerically the equation of motion between ri and r/. The 
temperature profile in the detached shell is constrained by 
the observed Hi line profiles. For the external part of the 
detached shell (between r/ and 12), Libert et al. (2007) as- 
sume a density profile that falls off as 1/r^ with the condition 
that the total mass (M_dt,cs) is given by the equivalent ISM 
mass enclosed in a sphere of radius r2 . Inside ri , Libert et 
al. assume a spherical wind in uniform expansion. 

Originally Libert et al. (2007) developed this model for 
Y CVn, a carbon star surrounded by a dust detached shell 
discovered by IRAS (Young et al. 1993a) and imaged by 
ISO (Izumiura et al. 1996). Hi line profiles obtained at dif- 
ferent positions on and around the central star were well 
fitted by such a model with a gas temperature decreasing 
from 1800 K, in ri, to 170 K, in r/. Libert et al. applied the 
same model to other sources, like RX Lep, an oxygen-rich 
star (Libert et al. 2008), and V1942 Sgr, another carbon star 
(Libert et al. 2010). In general, the gas detached shell model 
of Libert et al. (2007) accounts well for the H l emission of 
evolved stars, which are characterized by a narrow line pro- 
file (~ 2-4 kms~^) approximately centered on the stellar ra- 
dial velocity (Gerard & Le Bertre 2006). Indeed the emission 
from the detached shell where matter has been accumulating 
generally dominates the emission from the freely expanding 
wind region inside ri. The Hi emission from an unresolved 
spherical wind would be characterized by a rectangular line 




Figure 10. Schematic view of the a Ori detached shell. The 
termination shock is located in ri, the contact discontinuity in 
rj, and the bow shock in r2- CS stands for circumstellar material 
and EX, for external material. 



profile of width twice the expansion velocity (Gardan et al. 
2006). 

However, at high spatial resolution, the inner region 
(<ri) may be resolved. In such a case a spectrum obtained 
on the stellar position should show two horns separated by 
twice the expansion velocity. In addition, a central line, cor- 
responding to the detached shell should be observed in the 
middle, at the stellar radial velocity. Thus, for mass losing 
evolved stars, the Libert et al. (2007) model predicts a char- 
acteristic three-peaked profile at the central position, and 
a single-line profile away from it (>ri) with a maximum 
intensity close to the inner shell boundary. 

The merged VLA data show two peaks of emission cen- 
tered on a Ori at V*± Vcxp. Close to V*, one notes also a 
peak on the star (Fig. [T] upper panel) and a ring of emission 
around the star (Fig. [T] lower panel). Fig. [TT] shows spectra 
that have been extracted from the VLA merged data set 
at the position of the star and in a ring of internal radius, 
80", and external radius, 160". We expect that such a ring 
encloses most of the H l emission from the detached shell. 
In the upper panel of Fig. 1111 the predicted three-peaked 
profile is seen, with the peaks close to the expected veloci- 
ties. On the other hand a single peak at V* is visible in the 
spectrum shown in the lower panel. It thus seems that the 
VLA data contain information which, at least qualitatively, 
agrees with the Libert et al. (2007) model. However, the 
emission intrinsic to a Ori is partly hidden by the confusion 
due to the ISM emission, as evidenced by the negative peaks 
at 9kms~^. Similarly, the NRT data show an emission fea- 
ture centered close to V*. Due to its large beam, the NRT 
cannot detect the emission from the free-flowing wind of a 
Ori. 

In view of these results, we have run the Libert et al. 
(2007) model with parameters (distance, expansion veloc- 
ity, etc.) that apply to a Ori (see Sect. [2]). As the stel- 
lar effective temperature is larger than 2500 K, we assume 
that all hydrogen is in atomic form (Glassgold & Huggins 
1983). We assume also that 90% of the mass is in hydro- 
gen and 10% in helium, so that the mean molecular weight, 
^, is equal to 1.3. An expansion velocity of the free wind, 
vo = 14kms~^, implies a post-shock temperature in ri of al- 
most 5800 K (see Table [2]) . For the internal radius we select 



The detached shell of a Ori 9 



a value of 2' (=0.12 pc) in order to get a maximum Hi in- 
tensity at a position corresponding to the emission peaks 
identified in Fig.[71 The external radius was set arbitrarily to 
r2=3', and the ISM density to nn^lcm"^. The calculated 
Hi emission is convolved with the telescope beam profile. 
For the NRT, we adopt the response of a rectangular aper- 
ture of 160 mx 30 m (Gardan et al. 2006), and for the VLA 
a Gaussian response of FWHM = 34". 

We have tried several combinations of mass-loss rate 
and duration of the mass loss episode. The mass-loss rate 
is constrained by the intensity of the peaks observed at 
V*±Vcxp (i.e. -10 and -1-18 kms~^) on the star. The mass 
in atomic hydrogen in the detached shell (which can be es- 
timated by the emission feature observed at V*) constrains 
the product of the mass loss rate by the duration. We found 
satisfactory fits for values of M in the range 1.0 to 1.5x10"^^ 
M0yr~^ and values of the duration in the range 0.7 to 
1.0 xlO'^ years. 

In Figs.lHandllllwe present the results of our best com- 
promise that was obtained with M — 1.2x10^'' Mq yr~^ and 
a total duration of 0.8x10^ years. The parameters corre- 
sponding to this preferred model are summarized in Tabled 
A rapid cooling of the gas (with a temperature assumed to 
vary as a power law from 5800 K in ri to 220 K in Xf) is 
needed in order to reproduce the observed narrow line width 
(~ 3 kms~^) of the emission coming from the detached shell. 

In this spherical modeling we have adopted a systemic 
radial velocity of 3 kms~^, which gives the best agreement 
with NRT data. The H l emission from AGB stars is known 
to be sometimes shifted by ~l-3 kms~^ with respect to the 
CO emission which is believed to provide the most precise 
estimate of the stellar radial velocity (Gerard & Le Bertre 
2006). In general this shift is towards kms~^ (LSR), prob- 
ably an effect of the dragging of the circumstellar shell by 
the ambient ISM. In Fig. [TT] (upper panel) it is clear that 
the central peak is not centered at the average of the two 
horn velocities. Consistently the model predicts horns that 
are blueshifted by ~ 1 kms~^. 



6 DISCUSSION 

With the NRT and the VLA, we have detected an Hi source 
coincident with a Ori and relatively compact (~ 4' in right 
ascension). The emission is centered at Visr~ 3 kms~^, and 
has a narrow line profile (FWHM ~ 3 kms~^). In addition, 
the VLA channel map at 3.7 km s~^ shows a ring of emission 
encircling the position of the star and a tail extending south- 
west. Furthermore, the VLA data show emission lines at 
— 9kms~^ and at -1-18 kms^^, right at the position of the 
central star. Although our data are affected by confusion, 
these features can be clearly ascribed to the circumstellar 
shell expelled by the supergiant. 

These features can be explained qualitatively by the 
Libert et al. (2007) detached shell model. In particular, the 
spectrally narrow emission observed close to V*, which the 
VLA shows to peak along a ring encircling the position of a 
Ori at 3.7 kms~^, is characteristic of the emission expected 
from a detached shell. Using parameters which are represen- 
tative of a Ori, the model gives a fair agreement with the 
observations. The main result is that the detached gas shell 
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Figure 11. Upper panel: Hi profile obtained on a Ori with the 
VLA within an aperture of 34" diameter. Lower panel: Hi profile 
obtained by integrating the emission within a ring of internal 
radius, 80", and external radius, 160", and centered on the star. 
The dashed lines correspond to the modeled spectra discussed in 
Sect. [5] 



Table 2. Model parameters (d = 200 pc). The notations are as 
in Libert et al. (2007). In particular, tug is the formation time 
of the detached shell, y\-DT,CS is the mass of the circumstellar 
component of the detached shell, and Mjjy bx is the external 
mass accreted in the detached shell. 
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of a Ori can be accounted for by a mass loss episode of '^ 10^ 
years duration at a rate in the range l-2xl0~^ Mq yr~^. 

In Sect. 14.31 we noted that Bowers & Knapp (1987) de- 
tected only the horns at V*±Vcxp corresponding to the free 
flowing wind of a Ori, because they had selected large u-v 
spacings. They modeled these horns with a wind of 2.2xl0~^ 
Mq yr~^ (at 200 pc). In fact this agrees quite well with our 
model, because they adopted an envelope radius of 60" , half 
as much as our own ri estimate. The reprocessing of their 
full data set and the combination with D-configuration data 
revealed the stationary shell of size '-^4' that they missed. 

Here we want to stress again the importance of multi- 
scale observations, and in particular of the C-configuration 
data that we have retrieved from the VLA archive. Indeed, 
it is the improved angular resolution that allows us to char- 
acterize geometrically the detached shell. For three other 
sources that we have studied with the VLA in configuration 
D, Mira (Matthews et al. 2008), RS One (Matthews & Reid 
2007, Libert et al. 2010) and X Her (Matthews et al. 2011), 
we have not observed shells that are seen to be detached 
from the central stars. However the observed Hi line-profiles 
indicate that the stellar winds are slowed down by the exter- 
nal medium, and that the outer circumstellar shells should 
be quasi-stationary with an enhanced density. Thus we ex- 
pect that observations of these sources with a higher spatial 
resolution would reveal a detachment of the outer shells. 

The mass-loss rate that we find is smaller than the mass- 
loss rate used by Meynet & Maeder (2003) in the evolution 
of a 2O-M0 star which, among their models, gives the best 
agreement with the temperature and luminosity of a Ori. 
A better agreement on the temperature and mass-loss rate 
might be found for models with a lower initial mass of a Ori, 
but also with a lower luminosity, which would imply a dis- 
tance lower than 150 pc. We note however that, for the mass 
loss of cool stars, Meynet & Maeder (2003) use prescriptions 
from de Jager et al. (1988), which produce in their models 
smooth variations of the mass-loss rate with time. This is 
only indicative, as mass loss is known to be strongly vari- 
able in the progenitors of red supergiants, especially when 
they cross the instability strip. 

From observations, the mass-loss rate of a Ori is un- 
certain, and values ranging over a factor 10 have appeared 
in the literature. For instance De Beck et al. (2010), from 
a modeling of several CO rotational lines, derive a mass 
loss rate of ~ 5 x 10"'^ M0yr~^ (for a distance of 200 
pc), whereas Harper et al. (2001, 2009), from a model of 
the radio continuum emission, obtained 4.8 (±1.3) x 10~^ 
M0 yr~^ (for the same distance). This large range arises in 
part from tracers probing different regions of the circumstel- 
lar shell, and from a rate of mass loss that has likely been 
variable. Also the CO emission seems anomalously low in 
red supergiants which could lead to an underestimate of the 
mass- loss rate (Josselin et al. 1998). Interferometric data ob- 
tained at 11 /im by Danchi et al. (1994) seem to show the 
presence of several dust shells close (at 1 and 2") to the cen- 
tral star, suggesting a mass loss rate that is highly variable. 
We have adopted a value (1.2 x 10~® MQyr"^) which, in 
our model, represents an average over 8000 years. 

In our modeling the gas temperature decreases from 
5800 K (in ri) to 220 K (in r/) in a lapse of 72x10^ years. 
This translates to an average cooling time-scale of 2.2 x 10* 
years. The post-shock density is 56Hcm~^, and rises up 
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Figure 12. Simulation of the Cn line emission at various dis- 
tances from a Ori. For clarity, the spectra are successively shifted 
by 2K. 



to 1800 Hcm""^ at the contact discontinuity. Spitzer (1968) 
gives a cooling time of about 2.5x10 years for a gas of 
density lOHcm"^ and temperature 1000 K. Therefore the 
modeling is consistent with a quite reasonable hypothesis 
for the cooling. The cooling per unit mass is proportional 
to the density. From Dalgarno & McCray (1972) it should 
also depend on the electron density and on the tempera- 
ture. From our modelling we can estimate the temperature 
and the density as a function of the distance to the central 
star. The electron density is an unknown. However, hydro- 
gen should stay atomic within the detached shell, but species 
with a low ionization potential, such as carbon, are expected 
to be singly ionized (Huggins et al. 1994). 

The cooling lines from the detached shell are predicted 
to be centered at the same velocity as the Hi line (3kms~^), 
and to have a similarly narrow line-profile (~3kms~^). As 
an example we give the profile of the C 11 line at 158 ^m 
at various positions across the detached shell (Fig. I12|l . We 
assumed carbon to have an abundance of 2x 10"** relative to 
hydrogen and to be completely ionized beyond 6x10"'' pc 
(0.1') from the central star. The calculation has been done 
assuming a gaussian beam of 12" (FWHM), and adopting 
the same temperature profile as for hydrogen (Sect.[S]). This 
cooling line should be within reach of the IR Heterodyne 
Spectrometer on the SOFIA (Stratospheric Observatory for 
Infrared Astronomy). 

The efficient cooling that we have adopted implies that 
the detached gas shell in the model is geometrically narrow 
(r//ri ~ 1.2). This could be an explanation for the relative 
compactness of the Hi source (</; ~ 4') and for the squeezing 
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of the contours around the central star (except south-west) 
in the 2.4 and 3.7 kms~^ channel maps (Figs. [7] and |8]). 

The model that we are using is spherical and cannot ac- 
count for all the features observed in a Ori , which is known 
to be moving through the ISM at ~30kms~ towards the 
north-east. The plateau that is observed south-west in the 
3.7 kms~^ channel map could be the signature of a tail of 
material such as that observed behind Mira (Matthews et 
al. 2008) . This kind of feature is predicted by 3D hydrody- 
namic simulations of CSEs around mass losing stars moving 
through the ISM, such as those performed by Villaver et 
al. (2003) or by Wareing et al. (2007), and might hopefully 
contain information on past episodes of mass loss. 

The Hi emission that we observe in a Ori is concen- 
trated in a source of diameter ~ 4', and is thus found well 
inside the infrared emission reported by Noriega-Crespo et 
al. (1997) and by Ueta et al. (2008), and that Herschel has 
shown to arise from several thin shells of radii ~ 6' (Cox et 
al. 2012). These shells are clearly distinct from the classical 
detached shells of Young et al. (1993a). The mechanism re- 
sponsible for the far-infrared emission of these shells has not 
been identified with certainty. It could be emission by dust 
heated by collisions with the shocked gas and/or emission 
from forbidden lines of O I and C ll. However there is a con- 
sensus that this emission traces the bow shock(s) resulting 
from the motion of a Ori through the ISM. van Marie et 
al. (2011) have studied the coupling of dust and gas within 
a shocked wind with conditions typical of a Ori. They find 
that, apart for large grains (radius > 0.045 /xm), the dust 
remains bound to the gas in such a structure. 

We found however that the far-IR arc coincides with 
FUV emission and is associated with atomic hydrogen emis- 
sion (Figs. [S] and H|. The mass in atomic hydrogen com- 
pares fairly well with the estimate of the total mass of the 
far-IR arc by Noriega-Crespo et al. (1997), given the un- 
certainties in our Hi estimate and in the gas-to-dust ratio 
(100) adopted by Noriega-Crespo et al. Hi emission is found 
mainly at 7.6±1.3 kms~^, i.e. at the velocity at which inter- 
stellar emission peaks (see Fig. (J}. Thus the atomic gas that 
we detect in coincidence with the far-IR arc seems related 
to the surrounding ISM rather than to a Ori. Possibly, it is 
interstellar matter compressed in a precursor ahead of the 
bow shock. 

The ISM gas that crosses the bow shock at 
~30kms~^ should be heated at a temperature '^25 000K, 
and hydrogen may be ionized. Except to the south-west, we 
have found no unambiguous evidence of atomic hydrogen as- 
sociated with a Ori beyond 3' from the star. This suggests 
that hydrogen in the region between the interface and the 
bow shock remains ionized during the lapse of ~ 10^ years 
that we infer for the formation of the detached shell, or that 
it remains sufficiently warm that its emission is widened over 
a large range of velocities, preventing its detection. In these 
conditions, our Hi modeling cannot bring real constraints 
on the region of the detached shell which is filled by matter 
from the ISM. In particular, the value that we have adopted 
for the parameter r2 has no real significance. 

Conversely, the parameters that we find for the circum- 
stellar part of the detached shell could be used as constraints 
for the 3D hydrodynamic simulations of the Betelgeuse's 
bow shock that start to be developed (e.g. Mohamed et al. 
2011). In this specific case, we note that the interface is at 



^ 0.3 pc from the central star, which seems too large com- 
pared to our H I observations, that suggest a smaller value. 
Possibly the rate of mass loss adopted by Mohamed et al. 
(3 X 10"*^ M0 yr"^) is too large. The density of the ISM ma- 
terial through which a Ori is moving is also uncertain. 

Another possibility is that the far-IR arc seen at 6' from 
Q Ori is not the true bow shock related to the detached gas 
shell that is detected at 21 cm. The examination of the im- 
ages published by Cox et al. (2012) suggests that there are 
other structures inside the far-IR arc, in particular another 
arc at 3' north-east of a Ori. It is tempting to hypothe- 
size that this 3'-radius arc is the true bow shock related 
to the detached shell that we have modeled. In that case 
the 6'-radius arc which is also seen in H I and far-UV would 
be related to a more ancient episode of mass loss. Clearly 
spectral information on the far-IR emission of this arc and 
imaging in Hi at a higher spatial resolution are needed in 
order to elaborate on its true nature. The characteristics of 
the internal medium, between the detached shell and the 
far-IR arc, are also an issue. 



7 CONCLUSIONS 

Observations in the Hi line at 21 cm of a Ori have been 
performed with the NRT and the VLA. A source {(p ~ 4') 
has been detected at the position of the star. The emis- 
sion is dominated by a narrow line (FWHM ~ 3 kms~^) 
centered close to V*. The VLA data show that this source 
has the structure of a detached shell elongated towards the 
south-west. This shape agrees with the numerical simula- 
tions of a CSE moving through the ISM. Emission peaks at 
V*±Voxp are also observed by the VLA at the position of 
the central star. Presently the data bring no clear evidence 
of atomic hydrogen outside this source. 

We have developed a working model adapted to a 
Ori and its detached neutral gas shell that accounts for the 
spectral features that are observed and for the size of the 
source. In this model, the detached gas shell of a Ori results 
from the collision of a 1.2xl0~^ M0yr~^ stellar wind, ex- 
panding at 14kms~^, with the surrounding ISM for a lapse 
of ~ 10^ years. The gas is cooling from a post-termination- 
shock temperature of 5800 K to ~220 K, at the interface with 
ISM. 

The high level of confusion in the direction of a Ori 
forces us to be cautious. Nevertheless, by selecting different 
ranges of baselines, it is possible to obtain images revealing 
the detached shell of a Ori as well as H I emission coincident 
with the 6' far-IR/far-UV arc. Presently, we cannot conclude 
on the exact nature of the structure in arc discovered by 
IRAS. We suggest that the Hi emission might come from a 
precursor of the a Ori bow shock, or that the whole structure 
could be related to an ancient episode of mass loss. 
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